Quasiparticle Chirality in Epitaxial Graphene Probed at the Nanometer Scale 
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Graphene exhibits unconventional two-dimensional electronic properties resulting from the sym- 
metry of its quasiparticles, which leads to the concepts of pseudospin and electronic chirality. Here 
we report that scanning tunneling microscopy can be used to probe these unique symmetry proper- 
ties at the nanometer scale. They are reflected in the quantum interference pattern resulting from 
elastic scattering off impurities, and they can be directly read from its fast Fourier transform. Our 
data, complemented by theoretical calculations, demonstrate that the pseudospin and the electronic 
chirality in epitaxial graphene on SiC(0001) correspond to the ones predicted for ideal graphene. 



PACS numbers: 73.20.-r, 68.37.Ef, 72.10.Fk 

Graphene, one single layer of carbon atoms packed in 
a honeycomb structure, has two identical carbon atoms 
in each unit cell and thus two equivalent atom sublat- 
tices. This gives rise to an extra degree of freedom absent 
in conventional two-dimensional (2D) systems, and leads 
to exceptional electronic properties. The major conse- 
quence of the honeycomb structure is that the low en- 
ergy quasiparticles of graphene are described by a Dirac- 
like Hamiltonian. The energy spectrum is thus linear, 
and consists of two Dirac cones centered at the oppo- 
site corners K and K' of the Brillouin zone. This is, 
however, insufficient to understand properly the specific 
electronic properties of graphene. Due to the presence 
of the two C atoms per unit cell, the quasiparticles have 
to be described by two-component wave functions, each 
component specifying the (complex) amplitude on each 
atomic sublattice. This results in an additional degree 
of freedom known as pseudospin. The projection of the 
pseudospin on the direction of the momentum defines 
the chirality of the quasiparticules; in graphene the di- 
rection of the quasiparticle pseudospin is parallel (posi- 
tive chirality) or antiparallel (negative chirality) to their 
momentum [TJ [5]. Within one single Dirac cone, elec- 
trons (or equivalently holes) of opposite direction have 
opposite pseudospin, but the same chirality. In graphene 
bilayer with Bernal stacking, chiral Dirac fermions are 
also found, but they are massive and their chirality dif- 
fers from the one found in graphene monolayer [TJ [5] . 

In order to fully exploit graphene in future electronic 
devices it is of primary importance to assess the chiral- 
ity of the quasiparticles, since this feature is central to 
the transport properties. It is responsible for the new 
"chiral" quantum Hall effects (QHE) measured in ex- 
foliated monolayer and bilayer graphene [31 HJ [S]. In 
epitaxial graphene multilayers on the SiC(0001) face, 
Shubnikov-de-Haas oscillations and weak antilocalization 
effects consistent with the electronic chirality of "uncou- 



pled" graphene layers [51 [7] have been reported, but sur- 
prisingly the QHE has not been yet observed [5] . There- 
fore, a direct measure of the electronic chirality in a single 
epitaxial graphene layer is essential. Here we show that 
the quasiparticles pseudospin - and the associated chiral- 
ity - can be probed at the nanometer scale with a low 
temperature scanning tunneling microscope (STM). Our 
results on graphene monolayer on SiC(0001) are consis- 
tent with the electronic chirality of ideal graphene. 

Epitaxial graphene was grown on a 6H-SiC(0001) sub- 
strate by thermal desorption of silicon at high tempera- 
tures, in ultra high vacuum (UHV) jSJQTJ]. Parameters 
were chosen to obtain a mixed surface with a similar cov- 
erage of monolayer and bilayer (with Bernal stacking) 
terraces, ~ 50 nm wide. A carbon layer is present at 
the interface between graphene and the SiC. This layer 
is strongly bound to the substrate, and acts as a buffer 
layer with negligible interactions with the graphene layers 
[TTj . Although not metallic, the carbon buffer layer has 
electronic states that give a significant contribution to the 
STM images [UJJ [T^] . Our experiments were done using a 
home-made microscope operating at 4K in ultra-high vac- 
uum. High bias STM images, as shown in Fig. 1(a), allow 
discriminating monolayer from bilayer graphene terraces 
due to the different interface states contribution [TtJ] . 

Two groups have reported angle-resolved photoemis- 
sion spectroscopy (ARPES) measurements on epitaxial 
graphene layers on SiC(0001). They demonstrate that 
the graphene layers are electron doped. The measured 
difference between the Fermi energy (Ep) and the Dirac 
energy (Ed) is ~ 0.45 eV for a monolayer and ~ 0.30 
eV for a bilayer [T21 [21 [TS]. For graphene monolayer, 
a low energy linear dispersion has been found, a finger- 
print of Dirac-like fermions [T3J[T1]. However, different 
results arc obtained by the two groups concerning the 
crossing of the bands: An energy gap of 0.25 eV is found 
at the Dirac point by only one of the groups [2|- It has 



2 





\ K 


KV 








FIG. 1: (color online) (a) Constant current 90x43 nm 2 STM 
image of epitaxial graphene on SiC(OOOl), with two adjacent 
monolayer (ML) and bilayer (BL) graphene terraces. The 
spatial derivative of the image is shown, to enhance the cor- 
rugation due to the interface states which is higher on ML 
terraces [10]. Sample bias: +980 mV, tunneling current: 0.15 
nA. (b) Schematic Fermi surface for electron-doped ML and 
BL graphene. (c) Illustration of an intervalley scattering pro- 
cess, (d) 3D rendered 50x50 nm 2 constant current image of 
two adjacent ML and BL terraces, taken at low sample bias 
(+ lmV). A long wavelength scattering pattern is found on 
the BL terrace (left), and not on the ML terrace (right). Tun- 
neling current: 0.2 nA. (e) Schematic of forbidden intravalley 
backscattering for ML graphene. Red arrows sketch the pseu- 
dospin direction, (f) Schematic of intravalley backscattering 
for BL graphene. 



been proposed by the authors that the gap arises from 
a substrate-induced breaking of the sublattice symmetry 
of graphene. ARPES, also shows that the Fermi surfaces 
(FSs) of monolayer and bilayer graphene on SiC(0001) 
are essentially identical [131 U6j .Thev consist of two tiny 
circular pockets, of radius qp, surrounding the K and K' 
points [Fig. 1(b), six pockets are sketched to account for 
the symmetry of the Brillouin zone] . 

At low sample bias (a few mVs), the STM probes es- 
sentially the electronic states of the FS, and a constant 
current image reflects a map of the surface local den- 
sity of states (LDOS) at Ep [T7]. In a 2D system such 
as graphene, low bias images probe the quantum inter- 
ferences (QIs) arising from elastic scattering of surface 
quasiparticles at static defects [TH]. Such QIs are related 
to the Friedel oscillations in 2D metals. In a simple pic- 
ture, any impurity scattering between two states k and 



k ' of the FS should give rise to LDOS spatial modula- 
tions with wavevector k' — k , with a scattering proba- 
bility depending on the topology of the FS. For a circular 
FS of radius kp, scattering between states k p and states 
~ k p (termed backscattering in the following) are the 
most efficient processes due to enhanced phase space, re- 
sulting in LDOS modulations with wavector 2kp [181119] . 

For the FS sketched in Fig. 1(b), and depending on 
the impurity potentials present in the system, we may 
expect elastic electron scattering processes between the 
two non-equivalent Dirac cones at K and K' (i.e. in- 
tervalley scattering) as depicted in Fig. 1(c), together 
with elastic scattering processes between states located 
on the same Dirac cone (i.e. intravalley scattering). Fig- 
ure 1(d) is a constant current STM image taken at +2mV 
that includes two adjacent graphene monolayer and bi- 
layer terraces. A striking difference is measured in the 
LDOS patterns between the left and right terraces. A 
clear long-range modulation, of wavelength 5.2±0.3 nm, 
is found on the graphene bilayer terrace. This modula- 
tion has already been reported by Rutter et al. |12j . and 
is attributed to quantum interferences of wave vector 2qp 
due to intravalley backscattering between states ~~q p and 
states — ~q f within one single pocket [Fig. 1(f)]. Sur- 
prisingly, we do not find such LDOS modulation on the 
monolayer, which is puzzling if we keep in mind that the 
two systems have almost the same FS. 

The lack of LDOS modulation of wave vector 2qp for 
the monolayer implies that QIs associated with intraval- 
ley backscattering are virtually absent. It has been re- 
ported that wave function symmetries may lead to selec- 
tion rules that can prevent QIs, for instance in systems 
with strong spin-orbit coupling where the electronic spin 
governs the elastic scattering processes (2UJ, [2T] . In the 
present case, the absence of QIs in graphene monolayer 
is associated to the pseudospin instead of real spin. In- 
tuitively, this is because the pseudospin of the backscat- 
tered wave is opposite to that of the incident wave, so 
that their overlap should be zero [see Fig. 1(e)]. For 
a scattering potential with range larger than the C-C 
distance (which leads to scattering without pseudospin- 
flip), it has been shown that intravalley backscattering 
processes and the associated QIs are strictly forbidden in 
monolayer graphene [551 H3] ■ On the contrary, the elec- 
tronic chirality does not prevent intravalley backscatter- 
ing in bilayer graphene [Fig. 1(f)] [24] . 

For both monolayer and bilayer epitaxial graphene on 
SiC(0001), previous STM experiments have reported that 
atomic-size impurities give rise to short-wavelength mod- 
ulations of the LDOS, associated with intervalley scat- 
tering [lOj [12]. According to Ando et al. 22], such 
impurities, acting as short-range potentials, should also 
allow intravalley backscattering, which raises the ques- 
tion of possible QIs of wave vector 2qp. Recent theoret- 
ical studies show that for a general impurity potential, 
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FIG. 2: (color online) (a),(b) Low-bias STM images of 50 
nm wide monolayer (a) and bilayer (b) terraces. Sample bias 
and tunneling current are respectively +2 mV and 0.4 nA for 
(a), +4 mV and 0.13 nA for (b). (c),(d) Two-dimensional 
fast Fourier transform (FFT) maps of the STM images (a) 
and (b). (e) Central region of (c), showing no intravalley- 
backscattering related ring (the green arrow points out the 
position where such a ring should appear), (f) One of the 
outer pockets of (c). (g) Central region of (d), showing a 
clear ring-like feature of radius 2qF related to intravalley- 
backscattering. (h) One of the outer pockets of (d). Outer 
pockets shown in (f) and (h) are centered at the K (or K') 
point and result from intervalley scattering. 



LDOS modulations with wave vector 2qp should exist 
for graphene monolayer, but they should be strongly at- 
tenuated [251 HS1 121], with- a l/r 2 decay law instead of 
the l/r law found in conventional 2D systems. On the 
other hand, a l/r decay law is obtained for the QIs due 
to intravalley backscattering in a graphene bilayer, and 
also for the QIs due to intervalley scattering for both 
monolayer and bilayer graphene. In reference |25j one of 
us showed that the l/r 2 decay law should have a clear 
impact on the 2D Fourier transform (FT) of the LDOS 
modulations: while on bilayer graphene the efficiency of 
intravalley backscattering should give rise to a central 
ring of radius 2qp in the FT map similar to the one ob- 
served for conventional 2D systems, the strong attenua- 
tion of intravalley backscattering on monolayer graphene 
should generate a FT image with no central ring [25 . 

In the following, we present high resolution fast Fourier 
transform (FFT) of STM images recorded at low sample 
bias on 50 nm wide graphene terraces. We first show the 
real-space images, for graphene monolayer [Fig. 2(a)] 



and graphene bilayer [Fig. 2(b)]. Note that both im- 
ages exhibit a triangular pattern of periodicity ~ 1.9 nm 
which is related to the interface reconstruction [10] [12] , 
and which appears as a sextuplet of bright spots in the 
corresponding FFT images [Figs. 2(c) and 2(d)]. 

The central region of the FFT in Figs. 2(c) and 2(d) is 
related to intravalley scattering: a clear ring-like feature 
of average radius 1.2 nm" 1 is found for the bilayer [Figs. 
2(d) and 2(g)]. This radius value is in agreement with 
the value given in ref. [12 , and with the 2qp value de- 
rived from ARPES [13l[T4] . On the monolayer terrace, no 
central ring is found [Figs. 2(c) and 2(e)], despite the un- 
precedented momentum resolution obtained here. This 
has been checked on many different monolayer terraces 
(see also the EPAPS document [25J ), and we conclude 
that the result is robust and systematic. As discussed 
above, pseudospin and electronic chirality are the key 
parameters for the absence (presence) of the central ring 
in the FFT image of monolayer (bilayer) graphene. 

In the high frequency regions of the FFTs of Figs. 2(c) 
and 2(d), we find six outer pockets with ring-like shapes 
centered at K(K') points. They result from interval- 
ley scattering, associated to real space LDOS modula- 
tions with a (a/3 x \/3)i?30 o periodicity with respect to 
graphene [TUHH]. As shown in Figs. 2(f) and 2(h), the 
intensity of the high frequency rings in the FFT is not 
isotropic. The anisotropy is much more pronounced for 
graphene monolayer, for which a splitting of each ring 
into two arcs is found [Fig. 2(f))]. This splitting is dis- 
cussed in more detail in the EPAPS document [25] , 

The experimental results are compared with T-matrix 
calculation made along the lines of Ref. [23] (see also the 
EPAPS document [25]). Figure 3(a) is the FT of a cal- 
culated LDOS map at Ep for monolayer graphene with 
Ed = -Ef-0.45 eV, in the presence of a single localized 
impurity. It is clear that no central 2qp ring is generated 
in this case. This is the signature of LDOS modulations 
resulting from inefficient intravalley backscattering, and 
it is the most direct fingerprint of graphene's pseudospin 
and chirality |25j . This is in very good agreement with 
the experimental absence of the central ring [Fig. 2(e)]. 
Also, the high frequency split rings, resulting from inter- 
valley scattering, are nicely reproduced by the theoretical 
calculation. The splitting of the rings is another conse- 
quence of the quasiparticle chirality [28 . 

Figure 3(b) is an equivalent map as Fig. 3(a) but for a 
graphene bilayer, with Ep, = Ep-0.3 eV. From ARPES it 
is known that the two layers have a different electron den- 
sity, giving rise to a gap at Ep, of about 0.1 eV [15]. The 
asymmetry between the two layers has been included in 
the present calculation. Note the appearance of the inner 
ring of radius 2qp, as well as the presence of the outer 
rings which have not obvious splitting (a faint intensity 
remains perpendicularly to the TK and TK' directions). 
The agreement with the experimental FFT is excellent. 

Finally, we propose the use of STM to check any possi- 
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chirality of an ideal graphene sheet. Together with the 
previous evidence by ARPES of the low energy linear 
spectrum, a complete assessment of the Dirac nature of 
the quasiparticles in this system is obtained. 
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FIG. 3: (color online) (a-c) Theoretical calculation of the 
Fourier Transform (FT) of the LDOS at E F for different 
graphene layers with a single delta-function impurity, (a) 3D 
rendered FT map for graphene monolayer with a Dirac energy 
Ed = Ef— 0.45 eV. The calculation predicts no central ring of 
radius 2qp, showing the unefhciency of intravalley backscat- 
tering. (b) Same representation as (a) for a bilayer, showing 
a clear central 2qF ring. Ed = Ef— 0.3 eV. An asymmetry of 
~ O.leV between the two layers has been included, to account 
for the different electron density reported by ARPES (12,13). 
(c) Same representation as (a) for an asymmetric monolayer, 
with Ed — Ef— 0.45eV. A site energy difference of 0.25 eV is 
introduced between the carbon atoms of each unit cell, and 
the impurity is placed on a carbon atom with high site energy. 
The symmetry breaking between the two carbon sublattices 
restores the efficiency of intravalley backscattering as shown 
by the reappearance of a central ring in the FT map. 



ble breaking of graphene 's sublattice symmetry. To this 
end, we have calculated the FT of the LDOS map for 
a monolayer with the same parameters as in Fig. 3(a), 
but including an arbitrary sublattice asymmetry of 0.25 
eV, which is the gap value measured by Zhou et al. [T4"] . 
The calculation shows a clear inner ring in the FT, cor- 
responding to efficient QIs due to intravalley backscat- 
tering 28J . This dramatically illustrates the effect of the 
asymmetry on the FT STM data since this ring is ab- 
sent in the symmetrical case [Fig. 3(a)]. The data of 
Fig. 3(c) are for a defect on one sublattice, but we found 
that the shape of the ring depends on the location of the 
scattering potential [28] . Due to the large size of the im- 
ages it is likely that a distribution of defects is present 
in the experimental data. Calculations indicate that a 
large sublattice asymmetry (leading to a gap of 0.6 eV) 
should be readily detected in FT-STM images, but that 
for a moderate asymmetry (as for Fig. 3(c)) the inten- 
sity of the ring could depend on the defect configuration. 
Since we have taken images in many different surface ar- 
eas without detecting any ring, we believe that our data 
are rather consistent with a symmetrical monolayer, al- 
though a small asymmetry can not be totally ruled out. 

In summary, by analyzing the Fourier transform of low 
bias STM images, we have demonstrated that the pseu- 
dospin and the chirality of epitaxial graphene quasipar- 
ticles can be probed directly at the nanometer scale. We 
show that graphene monolayer on SiC(OOOl) exhibits the 
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